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A detailed analysis of the open charm effects on the decays of J/tp(ip') — > VP is presented, where 
V stands for light vector meson and P for light pseudoscalar meson. These are the channels that the 
_ so-called "12% rule" of perturbative QCD (pQCD) is obviously violated. Nevertheless, they are also 

the channels that violate the pQCD helicity selection rule (HSR) at leading order. In this work, we 
put constraints on the electromagnetic (EM) contribution, short-distance contribution from the cc 
annihilation at the wavefunction origin, and long-distance contribution from the open charm thresh- 
old effects on these two decays. We show that interferences among these amplitudes, in particular, 
«j ■ the destructive interferences between the short-distance and long-distance strong amplitudes play 



a key role to evade the HSR and cause the significant deviations from the pQCD expected "12% 
rule". 



(-| ■ PACS numbers: 13.25.Gv, 12.38.Lg, 12.40.Vv 
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INTRODUCTION 



Annihilation decays of heavy quarkonium have served as an important probe for the study of the 
perturbative QCD (pQCD) strong interactions in the literature [H-IH . In the bottomonium energy region, 
the non-relativistic approximation works well so that the annihilation of the bb can be regarded as a 
direct measurement of the properties of the bottomonium wavefunctions at the origin at leading order. 
For instance, for the 5* wave states, the annihilation matrix elements are proportional to the wavefunction 
at the origin, while for the P wave states to the first derivative at the origin. These simple relations have 
been broadly examined and found in good agreement with the experimental measurements in inclusive 
processes. They can be regarded as a direct test of the pQCD properties. Interestingly, although the 
mass of the charm quark cannot be regarded heavy enough, some of the leading pQCD relations are still 
well respected in inclusive transitions. A good example is the branching ratio fraction between ip' and 
J/yb: 



BR(ip' -> hadrons) ^ BR{iP' -> e+e") ^ 
= BR(J/ip hadrons) ~ BR(J/ip -> e+e~) ~ ' ' ( ' 

which is the so-called "12% rule" and the branching ratio fractions probe the ratio of the wavefunctions 
at their origins for the ground state J/tp and first radial excitation ip'. Note that in the above equation 
both branching ratios BR(J/ip hadrons) and BR(ip' — > hadrons) are referred to their light hadron 
decays. In fact, even for some of those exclusive decays, the above relation seems to hold approximately 
well. Such an observation, in contrast with the significant deviations in J/ip and ip' — > pir, has initiated 
tremendous interests in the study of transition mechanisms for J/tp and ip' — > pir, which is known as the 
so-called u pn puzzle". According to the Particle Data Group 2010 Q, the ratio for the pir channel is 
BR(ip' -> pn)/BR(J/ip -> pir) ~ (1.1 ~ 2.8) x 10~ 3 , which is much smaller than the pQCD expected 
value, i.e. ~ 12%. 

An alternative expression for the "p7r puzzle" is related to the power law suppression due to the pQCD 
helicity selection rule (HSR). As demonstrated in Refs. @, the decay of J/tp {ip') — > VP, where V 
and P stand for vector and pseudoscalar meson, respectively, should be strongly suppressed at leading 
twist. As a consequence, the branching ratio fraction is expected to be BR(ip' — > pn) / B R(J / i/j — > pir) ~ 
(Mj/^/M^i) 6 x (BR(ip' — > hadrons) / BR{J / ip — > hadrons)) ~ 6.2%, which is still much larger than the 
experimental observations. The significant violation of the pQCD HSR is nontrivial taking into account 
that quite many exclusive decay channels have approximately respected the 12% rule. 

Such a conflicting phenomenon has attracted a lot of attention from both experiment and theory in 
history. Even right now, the study of the "p7r puzzle" has been one of the most important physics 
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goals in the program of BESIII experiment [5|. In theory, this puzzle has also been broadly studied. 
Different explanations have been proposed in the literature, such as the color-octet model m, vector 
meson mixing @, HJ, final state interactions [1, admixtures of a vector glueball near 
intrinsic charm in light mesons [l3[ , light-quark mixing effects [TiJ , large survival decay of ip' via virtual 
charmonium state (l5j j , and interferences between the electromagnetic (EM) and strong interactions [l6l - 
[l9j . In the meantime, it has been realized that the "p7r puzzle" is not just restricted to the pn decay 
channel. It has also connections with the obvious charge asymmetries observed in ip' — > K*K + c.c. 
Therefore, it was conjectured that more general dynamic reasons should be investigated for J/ip (ip') — > 
VP [H,[2fl-Sl. 

It should be useful to recall the results of Ref. [3], where a global fit for J/ip(ip') — > VP is presented. 
The EM and strong transition amplitudes are parameterized out for all the decay channels, while among 
the strong transition amplitudes, the singly disconnected OZI (SOZI) processes and doubly disconnected 
OZI (DOZI) processes are further parameterized out. As shown in Ref. [ijj], there exists an overall 
suppression on the strong decay amplitudes of ip' — > VP, not just in the pir channel. Due to this 
suppression, the EM transition amplitudes become compatible with the strong decay amplitudes with 
which the interferences produce further deviations from the HSR-violating power law suppressions. This 
fitting result at least clarifies the following two issues: i) The same mechanism that suppresses ip' — > pn 
also plays a role in other ip' —> VP decays; ii) Such a mechanism does not affect much in J /ip — > VP as 
suggested by the charge asymmetries observed in K* K + c.c. These are important guidance for exploring 
mechanisms that would suppress the strong decay amplitudes in ip' — ¥ VP, but have less impact on the 
J/ip decays. 

During the past few years, we have been focussing on the study of mechanisms evading the HSR in 
charmonium decays. For charmonia below the open DD threshold, the HSR violating transitions are 
naturally correlated with the OZI-rule violations. As demonstrated in a series of studies [24- 28] , wc 
have shown that the intermediate D meson loops (IML) provide a natural mechanism for evading the 
OZI rule and hence the HSR in charmonium decays. The IML is introduced as a non-perturbative 
source of contributions. As iterated in Refs. [2(| HJ, [U, [25|, [29[, apart from the "p7r puzzle" the IML 
could be a key for understanding some of those long-standing questions in charmonium exclusive decays, 
e.g. the ip(3770) non-DD decay, large HSR-violating decay of r\ c — > VV, Ml transition problems with 
J/ip (i//) -» jr] c (r]' c ), etc. 

In this work, we provide a quantitative study of the role played by the long-distance IML in J/ip (ip') — > 
VP in association with the EM and short-distance SOZI transitions. Our purpose is to demonstrate 
that the IMLs as a non-perturbative transition mechanism are important for explaining the phenomena 
observed in J/ip (ip') — > VP, hence could be a natural solution for the long-standing "pir puzzle" and 
other puzzles in charmonium exclusive decays. As follows, the details of dealing with different transition 
amplitudes are given in Sec. II. The numerical results and detailed analysis are presented in Sec. Ill, 
and a summary in the last section. 

II. THE MODEL 

A unique feature with the VVP coupling is that at hadronic level the anti-symmetric tensor coupling 
is the only allowed Lorentz structure. Therefore, it can be understood that whatever the underlying 
mechanisms could be, they will contribute to the corrections to the anti-symmetric tensor coupling. 
Based on this, one can always make a general parametrization to the transition amplitude, 

Mtot = Mem + e i5 °(M sho rt + e ie M long ), (2) 

where Mem, M s hort and Mi ong are the amplitudes of the EM, strong short-distance and strong long- 
distance transitions. A phase angle 9 is introduced between the short and long-distance amplitudes, while 
the relative phase between the EM and short-distance amplitudes is 5q — 0° or 180°. It is reasonable 
to consider the trivial relative phase angles between the EM and short-distance amplitude. Meanwhile, 
the long-distance amplitude may carry a phase angle relative to the short-distance one due to hadronic 
wavefunction effects. Although the exclusive amplitudes for these three sources are obtained as real 
numbers, the relative phase angle 9 can lead to a complex coupling in J/ip(ijj') —> VP. We note that the 
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FIG. 1: The tree-level Feynman diagrams of EM transitions in J/ip(ip') — > VP. 



EM amplitudes for each decay modes carry intrinsic signs deduced in the quark model [16j . Our efforts 
as follows are to constrain these amplitudes and present an overall prescription for J/i/j(ip') — > VP. 



A. EM transition amplitudes 

The EM transition J/ip(ip') — > 7* — > VP turns out to be important in J/ip(tp') — > VP. In particular, it 
is the dominant contribution to those isospin-violating decay channels, i.e. J/ip(ip') — > pr), prj ', u>ir° and 
(f)ir° . This mechanism can be investigated in the vector meson dominance (VMD) model as presented in 
Refs. pHH. 

In Fig. Q] those three independent electromagnetic transition processes in the VMD are illustrated. 
The vertex couplings can be extracted from the experimental data for the decay widths of V — > jP (or 
P — > jV), and P — > 77. However, since the intermediate photon is off-shell, a form factor F(q 2 ) — 
^■em/(^em ~ I 2 ) 1S adopted for the EM transition amplitudes. The cut-off energy Aem is universal for 
both J/ip and ip' decays, and to be determined by experimental data for those isospin-violating decay 
channels. The EM amplitude can thus be expressed as 

M 

where p{k) is the four momentum of the initial vector charmonium (final light vector), and e(p) (e(fc)) is 
its corresponding polarization vector. In Tables HI \H\ and ITTll the EM vertex couplings are extracted with 
the up-to-date data from the PDG2010 



B. Short-distance transition amplitudes 



The short-distance contribution of strong interaction is mainly from the cc annihilation at the wavefunc- 
tion origin associated with hard gluon radiations. This is an SOZI transition and can be parameterized 
out in a similar way as in Ref. |18j |. We emphasize the cc annihilation at the wavefunction origin in this 
process. Thus, the HSR violation can be regarded as being produced by the non-negligible light quark 
masses in the hadronization process. The inclusive gluon annihilation part is thus guaranteed to scale 
with their lepton pair branching ratio fraction since the inclusive amplitudes would be controlled by the 
quarkonium wavefunction at the origin Q . This process distinguishes from the long-distance transitions 
via the IML where the cc annihilations would occur non-locally and probe the charmonium wavefunction 
away from the origin. Such a difference would allow us to treat the short and long-distance amplitudes 
individually and to avoid double-counting between these two mechanisms. A schematic diagram for the 
short-distance SOZI transitions is shown in Fig. [21(a). 
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TABLE I: The couplings g VlP = [12-kM^FCV -> 7-P)/|p 7 1 3 ] 1/2 or g Vl p = [4nM v F(P -> 7 V)/|p 7 | 3 ] 1/2 deter- 
mined by experimental data from PDG2010 Q]. 
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TABLE II: The couplings #p 77 = (32tiT(P -> 77)/M P ) 1/2 determined by experimental data from PDG2010 Q. 



9P~a 


Values 


r tot (keV) 


Branching ratios 


9^ii 


2.40 x 10~ a 


7.86 x 10~ a 


(98.823 ± 0.034)% 


9vTf 


9.68 x 10 _a 


1.3 


(39.31 ± 0.20)% 




2.13 x 10~' z 


194 


(2.22 ± 0.08)% 



The paramctrization of the short-distance amplitudes is outlined as follows [18j . First, the strength of 
the non-strange SOZI process is parameterized as 

= {(qq)v(qq)p\V \J/i>W)), (4) 

where Vq is the Zg decay potential of the charmonia into two non-strange qq pairs of vectors and pseu- 
doscalars via SOZI processes. But it should be noted that the subscript V and P here do not mean that 
the quark- antiquark pairs are the SU(3) flavor eigenstates of vector and pseudoscalar mesons. The am- 
plitude gj/^t^pi) is proportional to the charmonium wavefunctions at origin. Thus, it may have different 
values for J/ip and tp'. 

Considering the SU(3) flavor symmetry breaking, which distinguishes the s quark pair production from 
the u, d quarks in the hadronizations, we introduce the SU(3) flavor symmetry breaking parameter £, 

£ = ((qs) V (sq)p\Vo\J/^y))/gj^ m = {{sq) v {q- S ) P \V \J/^'))/g J/4 , m (5) 



TABLE III: The couplings e/fv = [3T v ^ e+e - /(2a e |p e |)] 1/2 determined by experimental data from PDG2010 0]. 



e/fv 


Values(xl0 _: ") 


rt ot (MeV) 


BR(V e + e") 


e/fp 


6.11 


149.1 


(4.72 ±0.05) x 10 _b 


e/U 


1.80 


8.49 


(7.28 ±0.14) x 10 _b 


e/U 


2.25 


4.26 


(2.954 ± 0.03) x 10" 4 


e / fj/i> 


2.71 


0.0929 


(5.94 ± 0.06)% 


e/U> 


1.62 


0.304 


(7.72 ± 0.17) x 10" J 
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FIG. 2: Schematic diagrams for the short-distance strong transitions in J/ip (tp') — > VP. Diagram (a) illustrates 
the SOZI process, while (b) is for the DOZI one. In both cases, c and c annihilate at the origin of the wavefunction. 



where £ = 1 is in the SU(3) flavour symmetry limit, while deviations from unity implies the SU(3) flavor 
symmetry breaking. In general, the value of parameter £ is around £ ~ fn/fic = 0.838, which provides 
a guidance for the SU(3) flavor symmetry breaking effects. For the production of two ss pairs via the 
SOZI potential, the recognition of the SU(3) flavor symmetry breaking in the transition is accordingly 

e = {{ss) v {ss) P \V Q \J/^'))/gj^ m . (6) 

For the J/ip and ip' decays into isoscalar final states, such as lot], wrf , 4>rj and cjn)', the DOZI transition 
as illustrated by Fig.^b) may also contribute. Although it is not apparent that the DOZI transition can 
be classified as a short-distance process, we can parameterized it out as follows, 

r = {{ss) v {qq) P \V Q \J/^{il}'))/gj/^) = ((qq)v(ss) P \V \ J/ip(ip'))/gj/^>) > ( 7 ) 

of which a small value \r\ << 1 would suggest a short-distance nature of this process. We mention that 
the DOZI process topologically does not double-count the long-distance IML transitions to be defined in 
the next Subsection. 

To take into account the size effects of the initial and final state mesons, a commonly adopted form 
factor is included, i.e. 

T(P) = \P\ l cxp(-P 2 /16/3 2 ) (8) 

where |P| is the three- vector momentum of the final-state mesons in the J/tp(ip') rest frame, and I is the 
final-state relative orbital angular momentum quantum number. We adopt j3 — 0.5 GeV, which is the 
same as Refs. (30l - [3^ |. At leading order the decays of J/ip(ip') — > VP are via P-wave, i.e. I = 1. 

The transition amplitudes for J/ip(ip') — > VP via the short-distance SOZI transitions can then be 
expressed as 

Ms(p°7T ) = M s {p + -K-)=Ms(p--K + ) = 9j/^ m T(V) 
M S (K*+K-) = M S (K*-K+) = M S (K*°K°) = M S (K*°K Q ) = gj Wr) ^(P) 
M s {uir,) = X r ,gj / ^, ) (l + 2r)F(P)+Y v V2Zrg J/iim F(P) 
M s (ujv') = X v ,g J/4 , m (l + 2r).F(P) + Y ri ,^rgj^ m F{V) 
M s {4>v) = X^rgj/^TiP) + ¥^,,,^{1 + r)£ 2 P(P) 

M S (H) = X^^rgj^T^+Y^gj/^il + r^HP), (9) 



G 



where X v (X v i ) and Y v {Y v i ) are mixing amplitudes between (uu + dd) / y/2 and ss components within the 
77 and 77' wavefunctions: 

77 = X v \uu + dd) / V2 + Y v \ss) , 

7/ = X v ,\uu + dd)/V2 + Y ri i\ss). (10) 

For the unitary 2x2 mixing, we have X v = — cos ap and X^ = —Y v = sinap with ap = 
dp + arctan(v / 2). The pseudoscalar mixing angle 9p is in a range of —22° ~ —13°. 

For the decays of J/ip(ip') — ► pn and K*K + c.c, the short-distance amplitudes are rather simple as 
listed above. For the decays into isoscalar final states, the situation would be complicated by the DOZI 
process and glueball mixing. There have been a lot of studies of the glueball mixing in the 77 and 77' 
wavefunction which can contribute to the isoscalar decay channels. However, in this analysis 

we do not consider the glueball mixing effects since the glueball components within 77 and 77' are rather 
small and need a delicate consideration. For the purpose of clarifying the role played by the short-distance 
and long-distance transition mechanisms in J /ip(ip') — > VP, we can leave the study of the glueball mixing 
effects to be considered in a differently motivated work [331 ] . 



C. Long-distance transition amplitudes via IML 



The IML transitions as a non-perturbative process seem to be a natural mechanism to evade the OZI 
rule and HSR in the charmonium decays [24|, HH, H3, [Hj]. The relevant effective Lagrangians for the 
charmonium couplings to the charmed mesons are as the following (36l . |37| : 

09 - 
C = i^Tr[R c - c H 2ll » d nH u ] + h.c, (11) 

where the S*-wave J/ij) and ip' charmonium states are expressed as 

Ro-c = (^) {Vl, - VcTo) (^) • (12) 
The charmed and anti-charmed meson triplet read 

Hu = (^)[X>r7,-2Vy6], (13) 
H 2l = [V*^-V a5 ](^iy (14) 

where T> and T>* denote the pseudoscalar and vector charmed meson fields respectively, i.e. = 
(U0W,D+(*),D+W). 

Consequently, the Lagrangian for the S'-wave J / ip and ip' is 

-ig^vv^V^^ - 9i n.T>er'#dr1> v <p£$prt + a vf), (15) 
The Lagrangians relevant to the light vector and pseudoscalar mesons are, 

C = -igvvv(p l d»V %] Vf -V«d*V ij &1) + ±g v . v . v e^ aP V*Sd v P i ^ a vf1 

+ igv*v*vVr^KK 3 (Vn) + 4»/d.„. v 23S(^V - VV»)p*J , (16) 
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with the convention e 0123 = 1, where V and V p are 3x3 matrices for the octet pseudoscalar and nonet 
vector mesons, respectively, 



V 



f^ + ^k * + K+ \ 




V K " -\/sV 

Based on the above Lagrangians, the explicit amplitudes in Fig. [3] can be obtained 

Mu p a \ t . P3XP3S \ 
vvv = -g^vvgwpe^pa^p^Q {-g\s-\ 2 — ) 

x {gvD*v9 u5 (P2 - Ps) ■ e - ±fv*v*vk 5 e v ) 
Mvftv* = -4g^v*vgv'V*vfv*V'v£p P aa£ aXKT eabcT^PiP%Pi\P3 K ,P%e b P3 

M Vtf)V = -g^VVgV'VVgVVV^ l ipaa^P P 1 P2 ( f{P2 ~ Pz) ' £ 

Mr u v a A ( i P2XP2S \ 
v , v , v = -g^VVgvVV3VVV^p,uafiP^ P2I {-g\5-\ 2 > 

m X>* 

x (4(P! - Ps)^ - (pi - pa) • cfl ff/5 + eJ( Pl - P2 ) 5 ) 
M v ,v*v* = gipvvgvvv£p,v aa pPiP'i{g-D>-D*vg pX (P2 -Pa) • e + ^fv-D*v^e x ) 

x (-.9A 5 + ^f^)«(pi-P2) 5 -e-(pi-p 2 )5^ + 4(Pi-P2r), (17) 

where p, fc, g are the four-vector momenta of the incoming charmonium, outgoing light vector, outgoing 
pseudoscalar, respectively, andpi, P2, P3 are the four- vector momenta of the intermediate charmed mesons 
as denoted in Fig. EJa). The subscriptions in the amplitudes denote the intermediate charmed mesons in 

the loops, and we have omitted the denominators, form factors, and integral measurement J f 2 ffl ^° keep 
the formulaes short. The following couplings are adopted in the numerical calculations [24l uEl [27l 128} : 



gibT>v = ^Wnhprn-D, g^vv — ~ , Sy>z>*z>* = g-ibw \ mv, Mx> = y/m-Dmv*, (18) 

Md V m f 

where g 2 = 2^"/^ ■ ana - m ^> ano - fi> ~ MeV are the mass and decay constant of J/tp. The relative 
coupling strength of -0' to J/ip, i.e. g^'vvl g.j/^vv — 1; i s included as a input. The light meson couplings 
to the charmed mesons are [38| 

2g g-D* vv 

gvvv — -r^/mvmv 1 gvv-v 



III r 



Po9v , fvvv ^gv m p 

gwv — gvw — — t=-, jvw = — —7=-, gv — -^i (.lyj 

V2 rnv V2 U 

where .9 = 0.59, /3 = 0.9, A = 0.56 GcV" 1 and / w = 132 MeV are adopted. 

The explicit amplitudes with different quantum number exchanges in the loops have been given in 
Eq. (fT7|) . For each decay mode the amplitude is dependent on the flavor component of the final state 
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light mesons. Thus, it is convenient to express the flavor-dependent amplitudes as 

1 

<2 ( 



M pri = Xn-dV^V } - [V+V+V-]) + c.c. 



M p , f = X n ^{[V V V°] - [V+V+V-]) + c 

M^o = ^{[V°V°V ]- [V+V+V-]) + c.c. 
M^o = 

M^a = ^{[V°V°V°] + [V+V+V-]]) + c.c. 

Mp+K- = [V°V°V+ +V-V+V ] 
1 

'2 
1 



M un = X v -{[V V°V a ] + [V+V+V-]) + c.c. 



= X V ^([V°V°V°] + [V+V+V-]) + c.c. 
M H = Y v [VfV;Vf]+c.c. 

Mw = Y v ,[vfv;vf]+c.c. 
M K , +K - = [v;vfv°] + [v v°vf] 

M K ,o RO = [v;vfv-] + [v+v-vt], (20) 

where X v (X v ') and Y v (Y v i) have been defined earlier, and the amplitudes of p-n+ , K*~ K+ and K*°K° 
have been implicated by their conjugation channels listed above. Note that the destructive sign between 
the charged and neutral meson loop amplitudes in those isospin- violating channels, such as prj, prj' , um a . 
The IML amplitudes for the 4>ir° channel vanish in the SU(3) symmetry limit. 

Since the IML integrals are ultra-violet divergent, an empirical tri-monopole form factor is introduced 

where to, is the mass of the exchanged particles and Pi is the corresponding four-vector momentum. As 
usual, Aj is parameterized into Aj = + chAqcd with Aqcd — 0.22 GeV denoting the typical low 
energy scale of QCD. 



III. MODEL RESULTS 
A. Analyzing scheme 



As mentioned earlier, all the underlying mechanisms in the VVP transitions would just contribute to 
the effective coupling constant. This feature, on the one hand, can provide advantages for disentangling 
different mechanisms, but on the other hand, may bring difficulties to the numerical fittings since the 
final results would only depend on the modulus of the summed amplitudes. Fortunately, the dynamic 
features of those different transition mechanisms as described in the previous section are useful for working 
out the parameter fitting scheme and disentangling the underlying mechanisms step by step. In Fig. 01 
we illustrate the relation between the EM and strong transition amplitudes (including the short and 
long-distance ones) by the addition of vectors in the complex plane. Our strategy of determining the 
amplitudes of those three transition processes is as follows: 

i) We treat the EM amplitude of each decay channel as a fixed vector in the complex plane pointing 
to the real axis as shown by Fig. [4] 

The EM amplitudes can be independently fixed by the data for those isospin-violating channels, i.e. 
J/ip(ip') —> prj, prf ', u!tt° and <jm°. The same parameter Aem — 0.542 GeV are then adopted for other 
decay channels as a reliable estimate of the EM amplitudes in the VMD model [l8j]. 
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J/W)(P) 






FIG. 3: Schematic diagrams for the long-distance IML transitions in J/ip (ij)') — > VP. In this case, the c and c 
annihilate by multi-soft-gluon radiations and can be described by intermediate charmed meson exchanges. The 
notations are T> = (D° , D + , Df) and T> — {D° ,D~ ,DJ). The light flavors of the charmed mesons are implicated 
by the final-state light mesons as shown in Eq. (|20[) . 



It should be noted that as discussed in Ref. [18| , the branching ratio fractions between ip' and J/ip decays 
into these isospin-violating channels are approximately within the range of 12% rule. This is an indication 
that for a single-mechanism-dominant process, the branching ratio fractions still serves as a probe for the 
wavefunction at the origin. In another word, if other mechanisms play a role, interferences among those 
processes would break down the pQCD relation. Such a scenario would be a natural explanation for the 
deviations observed in other channels. For instance, in the pir channel the interferences between the EM 
and strong amplitudes would lead to significant deviations from the 12% rule. Our focus in this work is 
to understand why the strong amplitude becomes compatible with the EM one in such an HSR-violating 
channel. 

ii) For the strong amplitudes including the short-distance and long-distance IML amplitudes, it is 
reasonable to impose that the EM and short-distance amplitude have a trivial relative phase since both 
probe the charmonium wavefunctions at the origin and both are real numbers. Moreover, for the short- 
distance SOZI amplitudes, we impose a constraint to require that their exclusive contributions should 
respect the 12% rule, i.e. the magnitude of the short-distance amplitudes can be treated as an input. 
Equation ([2]) can be rewritten as 



M t 



M 



EM 



6 ■A'l strong 5 



(22) 



where M. 



strong 



is the amplitude for the total strong transitions with a relative phase angle S relative to 



the EM one. By an overall fit of the experimental data [18j, the values of M. 



strong 



and S can be fixed for 



each channel. Then, with the fixed magnitude and direction for the EM and short-distance amplitudes, 
the decomposition of Eq. (|2"2"j) will allow us to determine the magnitude and direction of the long-distance 
amplitude as shown by Fig. 

We note that the overall fit of phase angle 6 in Ref. 
decays share the same value of S. Consequence of such an implication is a constraint on the magnitude 
and direction for the long-distance amplitudes in each VP channel. What we are going to examine in 
the following part is the range of the form factor parameter a in the meson loops, namely, whether all 



18] suggests that all the VP channels in J/ip or ip' 
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FIG. 4: Decomposition of the transition amplitudes for J/tp (if)') — > VP in a complex plane. The EM amplitude 
is assigned to point to the real axis while the short-distance amplitude carries a trivial sign difference to the EM 
one. The long-distance IML amplitude carries a phase due to hadronic effects. The final summed amplitude Mtot 
is to be compared with the experimental data. 

the VP channels share the same value of a at all. Confirmation of such a scenario should be evidence 
for the important contributions from the IML in J/tp and ip decays. 

Following the above procedure, we first consider pit and K* K + c.c. decay channels. Since the final 
state light mesons in these channels carry non-zero isospin, the short-distance transitions can only occur 
via the SOZI process while the DOZI process is forbidden. The analysis of these channels will then be 
able to expose the interfering feature of the IML. 

B. Parameters and results 

The parameters to appear in the analysis include the followings: i) the universal EM cut-off energy 
Abm = 0.542 GeV determined by the isospin- violating channels, ii) the short-distance transition strength 
9j/iji — 1-75 x 10~ 2 in the J/tp decay. This is an input for the short-distance amplitudes. It determines 
the short-distance transition strength g^> — 1.25 x 10~ 2 . Therefore, the exclusive contributions from the 
short-distance transitions still respect the 12% rule, iii) the form factor parameters aj/^ and a^i for the 
IML transitions which determine the long-distance coupling strengths, iv) the phase angles 6ji^ and 
between the EM and strong amplitudes in the J/tp and tp' decays, respectively, as defined in Eq. (|22|) . 

Other implicated parameters such as the SU(3) flavor symmetry breaking parameter £ = f-nj Jk — 
0.838, vertex coupling constants for the IML in Sec. Ill CI and flavor mixing angle 9p — —22° for r\ and 
rf , in principle, have been determined by independent processes. 

In Table IIV1 parameters adopted in the calculations are listed. As described in the above, some of 
those are treated as input, while the three parameters, i.e. the IML form factor parameter a, phase angle 
and SU(3) flavor symmetry parameter £, are fitted by the experimental data for J/tp and tp' — > pn 
and K*K + c.c, respectively. When the isospin zero decay channels are included, such as J/tp and 
tp' — s- uirj, uj-q 1 etc, the DOZI parameter r can also be fitted. We also list the x 2 f° r the J/tp and ip' 
decays, respectively. Briefly speaking, with the small number of parameters we can achieve a reasonable 
description of the overall experimental data, and the parameter uncertainties have been well constrained. 
The relatively large x 2 value for the J/tp decays are due to the relatively small experimental errors in this 
channel. One notices that the value of ocj/^ has relatively large uncertainties in comparison with that of 
the a^i . This is because the IML contributions are relatively small in the J/iJj decays, and insensitive to 
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TABLE IV: Parameters fitted respectively by experimental data for J/ip and ip' — > VP in our analysis scheme. 



parameter 








0.71 ±0.017 


0.92 ± 0.053 


r 


-0.201 ±0.006 


-0.097 ± 0.024 


e 


45.0° ±7.0° 


174.6° ± 2.9° 


a 


0.09 ± 0.03 


0.35 ±0.01 


x 2 


52.2 


9.17 



TABLE V: The effective couplings (in unit of GeV ) extracted from different transition mechanisms in 
J/ip (tp') — > pn and K* K + c.c. 



J/ip ->■ VP 


Qem 


Q short 




Qlong 




\Qstrong | 




\gtot\ 




p U 7T U 


— 2.26 x 10~ 4 


-2.36 x 10" 


-3 


-1.60 x 10" 


- b 


2.37 x 10" 


-a 


2.60 x 10 


-a 


K* + K~ 


-2.17 x 10~ 4 


-1.78 x 10" 


-3 


-1.45 x 10" 


-b 


1.79 x 10" 


-3 


2.0 x 10" 


3 




3.28 x 10~ 4 


-1.78 x 10" 


-3 


-1.45 x 10" 


- b 


1.79 x 10" 


-a 


1.46 x 10 


-3 


ip' -> VP 




-9.58 x 


-1.11 x 10" 


S 


-1.29 x 10" 


-3 


2.11 x 10" 


-4 


1.43 x 10 


-4 


K* + K~ 


-9.07 x 10~ s 


-1.07 x 10" 


-3 


-1.25 x 10" 


-3 


2.07 x 10" 


-4 


1.41 x 10 


-4 




1.37 x 10" 4 


-1.07 x 10" 


•A 


-1.28 x 10" 


-3 


2.33 x 10" 


-4 


3.58 x 10 


-4 



the form factor parameter. 

1. Isospin nonzero channels 

In Table [Vj the effective couplings extracted from different transition mechanisms in J/tjj and ip' — > pir 
and K* K+c.c. are listed. Comparing the coupling g s hort in these two channels, one notices that the SU(3) 
flavor symmetry breaking is at the scale of £ = 0.71 ~ 0.92 which is compatible with £ = f^/fx — 0.838. 
It shows that a small long-distance contribution from the IML will optimize the description of the data. 
The fitted form factor parameter aj/^p ~ 0.09 is adopted for all exclusive decay channels, which suggests a 
universal role played by the IML in J/ip VP. The small contribution from the IML is understandable 
since the mass of J/ip is much below the open charm threshold. Therefore, it docs not experience the 
long-distance IML effects in the transition. 

In J/ip — > VP, the relatively small EM amplitudes implies insignificant interferences between the EM 
and strong transition amplitudes. This feature is indicated by the relatively small charge asymmetries 
between J/ip -> K*+R- ± c.c. and J/ip -> K*°K° ± c.c. 

For ip' — > VP, the short-distance coupling strength is determined by the 12% rule relation. A proper 
description of the data leads to the determination of the long-distance IML amplitudes as listed in 
Table |Vj Since the mass of the ip' is much closer to the open charm threshold, the long-distance IML 
amplitudes become sizeable and play an important role in ip' — > VP. Given that the IML amplitudes are 
compatible with the short-distance ones in magnitude, the destructive interferences between the short 
and long-distance strong amplitudes have thus significantly lowered the strong transition amplitudes to 
be compatible with the EM ones. As a consequence, the further interferences with the EM amplitudes 
lead to the significant charge asymmetries between the branching ratios of ip' — > K*°K° + c.c. and 
iP' -> K*+K- + c.c. 

2. Isospin zero channels 

For the isospin zero decay channels, such as J/ip(ip') — > u)r\, ujrf ', cbru and cjyq', the DOZI transitions 
may contribute and two additional parameters have to be included [18j. One is the rj-rj' mixing angle 
ctp defined in Eq. (fTUj) and the other is the DOZI coupling strength r defined in Eq. We adopt the 
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commonly used value ap = 32.7° as an input, while treat r as a free parameter to be determined by 
the isospin zero decay channels. Meanwhile, all the other parameters determined in J/ip (-0') — > pn and 
K*K + c.c. are fixed. 

Eventually, it cannot be regarded as an overall fitting, and we do not expect a perfect description of 
the data for J/ip{ip') —¥ lot], lot]', <jyq, and cpr)'. This is mainly because the involvement of the DOZI 
mechanism and possible glueball mixing in the isospin zero channels should be considered in a more 
delicate way. Therefore, we only expect that those isospin zero channels to be described to the correct 
order of magnitude. 

In Table IVll the model calculations of the branching ratios of J/ip and ip' — >• VP are listed in comparison 
with the experimental values. The exclusive contributions from the EM, short-distance and long-distance 
IML transitions, and the combined strong contributions are also shown. The following points can be 
learned: 

i) For the isospin-violating channels, i.e. J/ip (ip') — > pr/, pr/ etc, the charged and neutral meson 
loops would cancel out exactly in the isospin symmetry limit. In other words, because of the isospin 
symmetry breaking, the mass difference between the u and d quark leads to ^ . As a result, 
the charged and neutral meson loops cannot completely cancel out, and the residue part will contribute 
to the isospin-violating amplitudes. Interestingly, we find that contributions from such a mechanism is 
much smaller than the EM transitions, which makes the isospin-violating channels ideal for the test of 
the 12% rule. In Table rvTfl the branching ratio fraction R is listed for all the VP channels. One can see 
that the 12% rule is reasonably respected in those isospin-violating channels. 

ii) For the final-state isospin nonzero channels, the systematic feature is that the long-distance ampli- 
tudes in J/ip — » pir and K*K + c.c. are negligibly small and the strong amplitudes are dominated by the 
short-distance ones. In contrast, the strong amplitudes are suppressed in ip' — > pir and K*K + c.c. because 
of the destructive cancellations between the short and long-distance amplitudes and become compatible 
with the EM ones. The observed branching ratio fractions are then further deviated from the 12% rule 
by the interferences between the EM and suppressed strong transition amplitudes in the ip' decays. We 
have shown that the suppression of the strong amplitudes is due to the open charm threshold effects 
via the IML transitions. The branching ratio fraction R is also listed in Table IVTll to compare with the 
data. Moreover, with the exclusive short-distance transition satisfying the 12% rule in the pn channel as 
a condition, other exclusive short-distance contributions also satisfy it fairly well. 

iii) For the final-state isospin zero channels, i.e. J/ip {ip 1 ) — > lutj and urj' etc, the experimental data have 
relatively large uncertainties and can be accounted for approximately to the same order of magnitude. 
Similar to the study of Ref. [l8j], the DOZI contributions are found necessary. In this analysis we do 
not consider the glueball mixing effects in the r\ and rj' wavefunctions since even though there might 
be glueball components within the rj and rj' , uncertainties caused by them may not be as large as other 
sources such as the DOZI contributions. In these processes, the EM amplitudes are relatively smaller than 
the strong ones. But the interferences among those strong transition amplitudes turn out to be sensitive. 
More delicate treatment for those isospin-zero channels are needed in further studies. One notices that 
the branching ratio fraction R can be reasonably accounted for except for the channels involving rj ' . This 
might be an indication that additional mechanisms should be considered. 

IV. SUMMARY 

By systematically analyzing the transition mechanisms for J/ip and ip' — > VP, we have shown that 
the long-distance IML transitions are crucial for our understanding of the long-standing "p7r puzzle" . 
Since the mass of ip 1 is close to the open charm threshold, its decays into VP are affected significantly 
via the IML transitions. In particular, the long-distance IML transitions provide a mechanism to evade 
the pQCD HSR and their destructive interferences with the short-distance amplitudes in the ip' decays 
cause apparent deviations from the "12% rule". The IML transition turns out to be a rather general 
nonperturbative mechanism in the charmonium energy region. Our analysis suggests that this mechanism 
should be present in all the decay modes. The same coincident cancelation between the short and long- 
distance amplitudes also causes large charge asymmetries between ip' — > K* + K~ + c.c. and K*°K° + c.c. 

It should be addressed that the open charm threshold effects via the IMLs can also contribute to the 
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TABLE VI: Theoretical results for the branching ratios of J/ip (ip 1 ) — > VP calculated in our model. The 
experimental data are from PDG2010 



BR{J/i: - 


> VP) 


EM 




short-distance long-distance 


strong 




total 
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X 


iu 


7 
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-3 
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-7 
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-3 
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X 
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3 
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pn 
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TABLE VII: The branching ratio fraction R = BR(i}}' -» VP)/BR( J /if) -> VP) given by our model (total) and 
the exclusive short-distance transitions. For the isospin-violating channels, the ratios are dominated by the EM 
transitions. The dash "-" means the transition is either absent or the data are not available. The experimental 
ratios are listed as a comparison. 



VP mode(%) 
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07T ,J 
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K*+R- +c.c. 


K* {> K { > + c.c. 
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11.98 


12.0 


20.92 


20.96 


18.41 


62.87 


23.13 


89.37 


total 


10.72 


10.80 
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10.27 


0.17 


0.18 


0.29 


3.45 


1.18 


20.89 


0.24 


64.65 


experiment 


7.40 


5.69 


3.0 






0.11 


0.19 


1.89 


< 0.63 


5.42 


2.41 


3.19 


- 16.50 


~ 41.38 


- 6.75 






~ 0.29 


~ 0.52 


~ 3.16 


~ 35.40 


~ 5.67 


~ 14.24 



process of ^(3770) — > VP. As shown in Refs. [ISIlSl, the IML mechanism can be a natural explanation 
for the sizeable ■0(3770) non-DZ) decay branching ratios observed in experiment (40T - |4^ ]. 

As a manifestation of the open charm threshold effects, the IML mechanism may also play an important 
role in other decay modes, such as J/ip(i/S) —> VS, VT, and PP etc. As the spin partners of J/i/jty'), 
the study of the ratio of BR{if c — » VV)/BR(r] c —> VV) should also be useful for clarifying the role 
played by the IML and provide some insights into the long-standing "p7r puzzle". Such a process has 
been investigated in Ref. [25| and recently updated in Ref. 26] . We expect that with the help of precise 
measurements of various decay modes at BESIII, the IML mechanism can be established as an important 
nonperturbative dynamics in the charmonium energy region. 
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